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| Abstract 

m ; 

For a multi-user interference channel with multi-antenna transmitters and single-antenna receivers, 
by restricting each transmitter to Gaussian input and each receiver to a single-user detector, computing 
the largest achievable rate region amounts to solving a family of non-convex optimization problems. 
Recognizing the intrinsic connection between the signal power at the intended receiver and the inter- 
ference power at the unintended receiver, the original family of non-convex optimization problems is 
ly-j ' converted into a new family of convex optimization problems. It is shown that, for such interference 

o ' 
in 
o 



channels with each receiver implementing single-user detection, transmitter beamforming can achieve all 
boundary points of the achievable rate region. 



■ Index terms — Gaussian interference channel, achievable rate region, beamforming 

O ' 

I. Introduction 

X 

The interference channel (IC) models a multi-user communication system in which each transmitter 
communicates to its intended receiver while generating interference to all unintended receivers. Deter- 
mining the capacity region of an IC remains an open problem except in the case of ICs with strong 
interference [1], [2]. To date, the best achievable rate region was established by Han and Kobayashi 
in [1], herein termed the HK region, which combines rate splitting at transmitters, joint decoding at 
receivers, and time sharing among codebooks. The HK region was simplified by Chong, Motani, Garg, 
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and El Gamal [3] and several computable subregions were also proposed in [4]-[6]. Etkin, Tse, and Wang 
[7, Theorem 1] proved that the HK region is within 1-bit of the capacity region of the Gaussian IC. The 
results in [6], [8] and [9], whose genie-aided approach is largely motivated by [7], established the sum- 
rate capacity of the two-user Gaussian IC in the low interference regime: when the interference power 
is below a certain threshold (referred to as noisy interference in [8]), the results assert the optimality of 
treating interference as noise at both receivers, i.e., each receiver should simply implement singe-user 
detection (SUD). In addition, even if the noisy interference condition is not satisfied, practical constraints 
often limit the receivers to implementing SUD. For example, the receivers may know only the channels 
associated with their own intended links. Under such scenarios, treating interference as noise at each 
receiver is more practical. 

There have been several recent studies of multiple-input multiple-output (MIMO) and multiple-input 
single-output (MISO) ICs [10]— [15]. The MISO IC describes, for example, the downlink communications 
of co-channel cells where the base stations have multiple antennas and the mobile stations have single 
antennas. The downlink beamforming problem has been well studied in [16]— [19]. Less well understood 
is the downlink transmission in the presence of interference, both in terms of its fundamental performance 
limits (i.e., capacity region), as well as in practically feasible transmission schemes. The assumption of 
multi-antenna transmitters and single-antenna receivers is motivated by the real world constraints where 
miniaturization of mobile units limits the number of antennas. In addition, the asymmetry in available 
resources at base and mobile stations favors systems where transmitters are tasked with heavy processing 
in exchange for reduced complexity at mobile units. Toward this end, we assume in the present work that 
each receiver implements SUD, i.e., it treats interference as channel noise. In a preliminary work [12], 
we showed that beamforming is optimal for the entire SUD rate region for a two-user real MISO IC. This 
result was used in [20] to characterize the beamforming vectors that achieve the boundary rate points 
on the SUD rate region. Later, the result in [12] was also used in [15] to derive the noisy-interference 
sum-rate capacity of the symmetric real MISO IC. In this paper, we generalize the result of [12] to 
complex multi-user MISO ICs. We note that the proof in [12] is applicable only to two-user real MISO 
ICs. 

There have been various studies concerning throughput optimization in a multi-user system under 
the assumption that each receiver treats interference as channel noise [21]-[26]. However, even for the 
simple scalar Gaussian IC, computing the largest achievable rate region with single-user detection at each 
receiver is in general an open problem [27]. Exhaustive search over the transmitter power is typically 
unavoidable due to the non-convexity of the problem. The difficulty is much more acute for the MISO IC 
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case as one needs to exhaust over all covariance matrices satisfying the power constraints. The complexity 
increases with the square of the number of transmit antennas, which renders the computation intractable. 
In this paper we propose an alternative way of deriving the optimal signaling for the SUD rate region for 
multi-user complex MISO ICs. Our approach is to convert a family of non-convex optimization problems 
for the original formulation to an equivalent family of convex optimization problems. What is more 
significant is that, given that each receiver implements SUD, all boundary points of the rate region can 
be achieved by transmitter beamforming. 

The rest of the paper is organized as follows. In Section II, we use the two-user complex MISO IC as 
an example to explain the basic idea of problem reformulation. We show that beamforming is optimal 
for the SUD rate region of such channel. The closed-form rate region is also presented. These results are 
generalized to multi-user complex MISO ICs in Section III. We prove that beamforming is also optimal 
for m-user complex MISO ICs with m > 2. Based on this result, we use the three-user MISO IC as an 
example to show how to obtain the SUD rate region for an m-user MISO IC. Numerical examples are 
provided in Section IV. We conclude in Section V. 

Before proceeding, we introduce the following notation. 

• Bold fact letters, e.g. x and X, denote vectors and matrices respectively. 

• (-) T and (-)t denote respectively the transpose and the Hermitian (conjugate transpose) of a matrix 
or a vector. Consequently the Hermitian of a scalar is its conjugate. 

• I is an identity matrix, is an all-zero vector or matrix depending on the context, and diag(x) is a 
diagonal matrix with its diagonal entries the same as that of the vector x. 

• X y means that X is a positive semi-definite Hermitian matrix. 

• tr(X) and rank(X) denote the trace and the rank, respectively, of the matrix X. 

• (x)i denotes the ith entry of vector x, (X)^ denotes the ith row and jth column entry of matrix 
X, and X mxn means that X is an m x n matrix. 

• \x\ is the absolute value of a scalar x, and ||x|| is the norm of a vector x, i.e., ||x|| = Vx^x. 

• Z (x,y) denotes the angle between two real vectors x and y, and Z (x,y) G [0, ir}. If both x and y 
are non-zero, then Z (x,y) = cos -1 pjp^j] ■ Otherwise we let Z (x,y) = § for convenience. 

• E[] denotes expectation. 

• sign(x) is the sign of a real scalar x, i.e., 

[-1, if x<0, 
sign(x) = < 

1, if i>0. 
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• Re(-) denotes the real part of it complex argument. 

II. Two-user MISO IC with single user detector 
The two-user Gaussian MISO IC is illustrated in Fig. 1 and the received signals are defined as 




Fig. 1. The two-user Gaussian MISO IC. 



Yi =h\xi+ h\x 2 + Ni 
and Y 2 = h\xi + h\x 2 + N 2 , 



(1) 



where X\ and x 2 are transmitted signal vectors of user 1 and user 2 with dimensions t\ and t 2 , respectively; 
Y\ and Y 2 are two scalar received signals; N\ and N 2 are unit variance circularly symmetric complex 
Gaussian noises; h\ and hz are t\ x 1 complex channel vectors; and h 2 and are t 2 x 1 complex 
channel vectors. The power constraints at the transmitters are respectively tr(Si) < Pi and tr(S2) < P 2 , 



where S\ = E 



x\x\ 



E 



x 2 x\ 



. We assume that the transmitted signals x\ and x 2 are zero-mean 
Gaussian vectors. Each transmitter knows all the channel vectors. Each receiver knows only the channel 
vector from its transmitter: receiver 1 (resp. 2) only knows (resp. h±). Each receiver decodes its own 
signal while treating the interference from the other user as noise. The boundary points of the achievable 
rate region for this channel are characterized by the following family of optimization problems: 



max (i±R± + fJ- 2 R 2 
subject to Ri = log ^1 + 

R 2 = log ^1 + 

tr(S0 < Pi, 
Si^O, i 



h\S\hi 
l + h\S 2 h 2i 

h\s 2 h,4 
l + hls 1 h 3/ 



1,2, 



(2) 
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where < fi\ < oo and < fi2 < oo. We define the SUD rate region of a two-user MISO IC as 

[J {R± and R2 optimal for problem (2)} . (3) 

all /ii,/i2>0 

Problem (2) is a non-convex optimization problem. For each \i\ and ^2 pair, all possible Si and S2 
must be exhausted over to find the solution of problem (2). To obtain the entire SUD rate region, one 
has to go through this exhaustive search for all the \i\ and [12 pairs. This exhaustive search is prohibitive 
when t\ and £2 become large. 

In the following, we first convert problem (2) into a family of convex optimization problems, and then 
obtain their closed-form solutions. 

A. Problem reformulation 

We first define the following optimization problems: 

max /i[Si/ii 
subject to /I3S1/13 = z\ 

tr(Si)<Pi, SibO, (4) 

and 

max /»|S2^4 
subject to /4S2/12 = z\ 

tr(S 2 )<P 2 , S 2 ^0. (5) 
In order for problems (4) and (5) to be feasible, we require 



0<^i< max x/htS^ = v/PiHM (6) 

It S; : P 1 ,S -0 V 



and < Z2 < max a/zJ,^^ = V^ll^ll- (7) 



tr(S 2 )<P 2 ,S 2 hO v 2 " 

We now establish the equivalence between problem (2) and the above two optimization problems. 
Lemma 1: For any non-negative scalars \i\ and \i2, the optimal solution and S2 for problem (2) is 
also an optimal solution for problems (4) and (5) with z\ = z* 2 = h\S\h^ and z\ = z* 2 2 = hls^. 
Proof: Problem (2) is equivalent to the following optimization problem for the same \i\ and \i2. 



max 



A*l I 1 + /, ...2 + log 1 + -T—2 



1 + Z 2 J V 1 + Z V 

subject to ^gSi^3 = z^ 2 , h\S2h2 = z 2 2 , 

tr(S;)<P, SihO, * = 1, 2. (8) 
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The equivalence is established as follows. First, the maximum of problem (2) is no smaller than that of 
problem (8), since problem (8) has extra constraints h\Sihs = z\ 2 and h\$>2h2 = z| 2 . On the other hand, 
the maximum of problem (2) is no greater than that of problem (8), since S\ and S2 are also feasible 
for problem (8), which are the optimal solutions for problem (2). Therefore, problems (2) and (8) are 
equivalent. We now recognize that problem (8) is equivalent to problems (4) and (5) with z\ = z\ and 
%2 = %2> which can be solved individually. ■ 
We remark that the optimization problem (8) can not be solved independently as the constraint 
parameters z\ and Z2 depend on the unknown optimal covariances. That is, unless the optimal SJ and S?; 
of problem (2) are obtained, the equivalent optimization problems in the form of (4) and (5) cannot be 
parameterized. However, this problem reformulation becomes especially powerful when we need to find 
the entire achievable rate region (or its boundary points) and to study the optimal signaling structure. 
Even though one cannot solve any individual optimization problem (2) by the corresponding problem 
(8), Lemma 1 establishes the following crucial fact that enables us to obtain the entire SUD rate region 
without explicitly solving (8): 



where the left-hand side denotes the collection of all the optimal solutions of problem (2) found by 
exhausting over m and fi2, and the right-hand side denotes the collection of all the optimal solutions of 
problems (4) and (5) found by exhausting over z\ and z 2 . Since the SUD rate region is determined by 
the left-hand side of (9), Lemma 1 successfully converts a family of non-convex optimization problems 
(2) into a family of equivalent convex optimization problems (4) and (5). 

To be more precise, instead of solving the family of non-convex optimization problems by exhausting 
over hi and ^2, one can instead solve the family of convex optimization problems by exhausting z\ and 
Z2 over the range specified by (6) and (7). We now proceed to obtain closed-form solutions for problems 



B. Optimal Solution 

By symmetry, we need only solve problem (4). Assume the singular-value decomposition (SVD) of 
hi is 




(9) 



(4) and (5). 



hi = 



Ui 



i = 1,...,4, 



(10) 
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/ill 

/3 



where /in = and j3 is a (ii — 1) x 1 vector. We have the following lemma. 

Lemma 2: Assuming the optimization problem (4) is feasible, the following SJ's are 
• If h\ and h% are linearly independent (consequently \\hs\\ / 0, and ||/3|| / 0), then 

Si=7i7 f i, 



and the achieved maximum is 



z\ 



h\h± 



\h,\\ 2 



+ 



\hi\ 



h^hi 



Pi 



where 



7i = u 3 



Zl 



h, 



k* Pi 



r l\ 



and k 



11 



hu 
1, 



z 2 

IIM 2 ' II/3II 
if /in / 

otherwise. 



• If h\ and /13 are linearly dependent and \\h^\\ / 0, then 



St = 



7 2 
Z l 



\h 3 



and h[sihi 



\h,P 



12^3^3 



• If ||/i3|| = (hence z\ = 0), then 



St 



Pi 



\\hi\\ 2 

and /i^S^i = Pi||/m || 2 



/ii/i{, 



Moreover, for all above cases, we have 



\h,P 



rank(S^) < 1. 



The proof is given in Appendix A. 



July 3, 2009 



8 



Lemma 2 shows that, for a fixed interference power z\ , transmitter beamforming maximizes the received 
signal power 1 . If hi and hs are linearly independent, the quadratic constraint of (4) defines a set of 
beamforming vectors whose projections on h 3 have equal length z\. Among all these vectors, the one 
that has the largest length of the projection on h\ is the optimal beamforming vector. 



C. The SUD rate region of a two-user MISO IC 

With Lemmas 1 and 2, we obtain the SUD rate region of a two-user MISO IC. 
Theorem 1: The SUD rate region of a two-user MISO IC with complex channels is 

Pi||/ii|| 2 sin 2 (^+Vi; 



U 

V>ie[o,-f-e 
</>2e[o,f-0 2 ] v 



Ri < log (l + 

R 2 < log ( 1 + 



where 



cos 



h^hi 



and #9 = < 



l*l|| " 11^3 



COS 







\\h2\\ ■ ||M 



l + P 2 \\h 2 f sin 2 i) 2 
Pall^fsin^^ + ifo) 
l + Pi||/i 3 || 2 sin 2 Y>i 



if||fei||^0, ll^sll^O 
otherwise, 

if||h 2 ||^0, ||fc4||^0 
otherwise. 



(21) 



Furthermore, the boundary points of the rate region can be achieved by restricting each transmitter to 
implement beamforming. 

Proof: We first assume that h\ is linearly independent of h.3, and that h 2 is linearly independent of h^. 
Define 



<j)\(z\ ) = tan 1 



Pi 11*3 



v2 ' 



then (13) becomes 

ftlSjfci = Pi ll^i || 2 sin 2 {6 1 + Mzi)) ■ 
Similarly, the maximum of problem (5) is 

h\S* 2 h 4 = P 2 \\h4 2 sin 2 (6 2 + <f> 2 (z 2 )) , 

'in Appendix A, we show that when h\hs = 0, there exist some matrices that are not beamforming matrices but still maximize 
problem (4). This also happens when hi and hi are linearly dependent. However, these does not contradict our conclusion that 
beamforming is optimal. 
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where 

Therefore, the achievable rate region determined by problem (2) is 

' Ri< ^flWWftiM,)]' 

V 1 + ^2 

P 2 < log ( H r~rz2 



u 



zfe[o,Pi\\h 3 \\ 
zle [o,P 2 \\h 2 1| 



(22) 



On defining fa = tan -1 Y ^pfft^ ? an( * ^ 2 = tan 1 \j p^h?\\^ ~ ' ^ e mter f erence P ower an d the 
useful signal power caused by transmitter 1 are given respectively by 

z\ = Pi ||^3 1| 2 sin 2 ^, (23) 

and 

fcJSjfci = Pi||/ii|| 2 sm 2 (#i + ^i). (24) 

Similarly, the interference power and the useful signal power caused by transmitter 2 are given respectively 
by 

4 = -P 2 ||/i2|| 2 sm 2 V 2 , (25) 

and 

h\S* 2 h 4 = P 2 \\h 4 f sm 2 (9 2 + fa). (26) 

Since fa varies continuously in [0, §] as zi varies in [0, P||/&3|| 2 ] (and similar for fa), the region in 
(22) is the same as 

Pi||fci|| 2 sin 2 (9 l + fa) 




Pi < log 1+ , . , 

V 1 + P 2 ||^2f 5111^2 / V ni , 
P 2 < log 1 + ~ pill. ||2 ' 2 / 

V 1 + Pi||/*3|| 2 sm^i 

When G [| — 6i, ^] , the useful signal power Pi \\h\ || 2 sm 2 (#i+V>i) or P2II/14II 2 sin 2 (^2+V'2) decreases 
as the interference power Pi || J13 1| 2 sin 2 fa or P2 \\h 2 1| 2 sin 2 fa increases. As such, the rate pairs associated 
with fa E [f — Oi-, f ] are interior points of the set (27). Therefore, (27) can be simplified into (21). 

In the cases where h\ and h 3 are linearly dependent, (21) is still the SUD rate region. This is due to 
the fact that (17) can be expressed in (24) when 6\ = 0; and (19) can also be expressed in (24) when 
6>i = and fa = f . 
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For each choice of ipi and fa, the corresponding S* can be obtained from (12), (16) or (18), and 



The rate region in (21) is characterized by ipi and ip2 which also determine the interference powers 
(23) and (25) at the two receivers. Compared to the original problem (2) which is characterized by the 
slopes (— j^) of the boundary points and requires the solution of a family of non-convex optimizations, 
Theorem 1 gives closed-form solutions by problem reformulation. Moreover, Theorem 1 shows that to 
achieve the rate pairs on the boundary of the SUD rate region, the transmitter can restrict itself to a 
simple beamforming strategy. 

In the SUD rate region (21), we point out several special rate pairs. 

• tpi = 4 , 2 = 0. This corresponds to a zero-forcing (ZF) beamforming rate pair as both transmitters 
generate no interference to the unintended receiver. Thus R x = log(l + Pi||/ii|| 2 sm 2 6>i) and R 2 = 
log (l + P2||^4|| 2 sin 2 82) are the maximum ZF beamforming rates. In general, this rate pair is in 
the interior of the SUD rate region. 

• tpi = 0, tp2 = § — 02- This case shows that user 1 can communicate at a rate no greater than 
Ri = log (l + Pi||/ii|| 2 sin 2 8\/ (l + P2||^2|| 2 cos 2 8 2 )) when user 2 is at the maximum rate R 2 = 
log (l + P2H/14H 2 ). This corresponds to a corner point on the rate region. 

• ipi = § — 8\, ip2 = 0. This is the other corner point of the rate region. User 2 can communicate at 
a rate no greater than R 2 = log (l + ^2 1 1 1 1 2 sin 2 8 2 / (l + Pi||/i3|| 2 cos 2 #1)) when user 1 is at the 
maximum rate R\ = log (l + Pi||/ii|| 2 ). 

When h\ and /13, and /12 and are respectively linearly independent, both transmitters use all their 
power to achieve the largest rates. However, if either of the above two vector pairs is linearly dependent, 
the transmitters do not necessarily use all the power. An example is the scalar Gaussian IC. 

Lemma 3: [28, Theorem 6] If h\ = = 1, h,2 = \fa and /13 = y/b, then the maximum SUD sum 
rate is 



Therefore, the maximum SUD sum rate for a scalar IC is achieved by letting both users use all the 
power, or letting one user use all the power while keeping the other user silent. The contrast between a 
scalar IC and a MISO IC with linearly independent channels is largely due to the existence of and the 



similarly for S|. 



R 3 = max {/(Pi, P 2 ), /(0, P 2 ), /(Pi, 0)} 



(28) 



where 
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interplay between the spatial diversity and multi-user diversity of a MISO IC. There is a tradeoff between 
the power of the intended signal at its own receiver and the interference power at the other receiver. For 
a scalar Gaussian IC, these two signals overlap in the same subspace and are proportional to each other, 
i.e., the channels are always linearly dependent. An increase of the intended signal power always results 
in an increase of the interference power (see (23) and (24) when 6\ is 0). The optimal tradeoff is achieved 
by choosing the appropriate power at the transmitters. As shown in (28), the optimal tradeoff does not 
necessarily require that both users use all the power. But for the MISO IC with linearly independent 
channels, the intended link and the interference link are in non-overlapping subspaces. Therefore, the 
optimal tradeoff is achieved by choosing the optimal beamforming subspaces while using all the power. 

Fig. 2 is an illustration of the beamforming vector of MISO IC. For simplicity, the channel vectors 
hi and h3 are assumed to be 2 x 1 real vectors with unit lengths. The angle between hi and Ii2 is 6. 
The disc with radius y/P contains all possible beamforming vectors that satisfy the power constraint. E 
and F are on the circle, the projections of vectors OE and OF on I13 both have length z. Then all the 
vectors on the line segment EF satisfy the power constraint P and the interference constraint z. Among 
those vectors, OE has the greatest length of projection on hi. Therefore, OE is the optimal beamforming 
vector 7 given the interference constraint z. It can be shown that the angle between OE and EF is tp 
and the length of the projection of OE on hi is y/Psm(6 + tp). 

The reduction of the non-convex optimization problem (2) to the equivalent optimization problem (4) 
is obtained by fixing the interference power while maximizing the useful signal power. This method 
is equivalent to fixing the useful signal power while minimizing the interference power. This requires 
solving the following optimization problems: 

min /I3S1/13 

subject to h\S\hi = zf , 

tr(Si)<Pi, Si^O, 

and 

min ^^2^2 
subject to h\S2h,4 = z\. 

tr(S 2 )<P 2 , S 2 ^0. 

The above two problems can be solved in the same way as problem (4), and the rate region can be similarly 
obtained. For these two proposed methods, the constraints are imposed either on the interference powers 
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Fig. 2. A geometric explanation of beamforming in the MISO IC. 

or the useful signal powers. One can also combine these two methods. For example, we can impose 
a constraint on the interference power caused by transmitter 1 while maximizing the useful power for 
receiver 1, and in the meantime, impose a constraint on the signal power on receiver 2 while minimizing 
the interference power caused by transmitter 2. 



D. Interference-limited SUD rate region 

As Theorem 1 is obtained by examining the relationship between the interference power and the useful 
signal power, we can easily apply it to the MISO IC under interference power constraints. 

Theorem 2: For the MISO IC defined in (1) with \\h2W 7^ 0, \\h3W 7^ 0, and with two additional 
constraints < Q\ and ^^2^2 < Q2 on the interference powers, the SUD rate region is 

Pi||/H|| 2 sin 2 (0i+V>iy 




l + P 2 \\h 2 \\ 2 sin 2 V2 
P2||M 2 g in ( g 2 + ^2)' 
1 + Pi||ft 3 || 2 sin 2 V>i 



(29) 
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where 




sin 1 



Qi 



aq ■ J n a . -i Q 

and 62 = mm < — — 02, sm 



L 2 Y P2\\h 2 \\ 2 j ' 

The proof is straightforward from Theorem 1. 

When Qi > .Pi ||^. 3 || 2 cos 2 6*i and Q 2 > J P2||^2|| 2 cos 2 6» 2 , (29) is exactly (21). Therefore, when the 
interference constraints are larger than the above thresholds, these constraints do not change the SUD 
rate region. Another extreme case is that in which neither user is allowed to generate interference to the 
other user. This is also the ZF rate region which is a rectangle determined by Qi = Q 2 = 0. 

III. Multi-user MISO IC with single-user detection 

In this section, we generalize our study of the two-user case to the general multi-user MISO IC. The 
key is, again, the problem reformulation as illustrated in Lemma 1 for the two-user case. For the general 
m-user MISO IC, we prove the optimality of beamforming with an SUD receiver. We then given the 
explicit description of the SUD rate region for a three-user MISO IC, and generalize it to the m-user 
case. 

A. Optimality of beamforming for an m-user MISO IC 
Define the received signal of the zfh user, i = 1, ■ ■ ■ , m, as 

m 

Y i = Y J h) i x 3 + N t , (30) 

where Xj is the tj x 1 transmitted signal vector of user j, hji is the tj x 1 complex channel vector from 
the jth transmitter to the ith receiver, and Ni is unit variance circularly symmetric complex Gaussian 



. As with the two-user case, 



noise. The power constraint for user i is tr (Sj) < Pi where Sj = E 
the input signals xfs are all zero-mean Gaussian vectors and each receiver treats interference as noise. 
Lemma 1 can be easily extended to multi-user MISO ICs, as follows. 

Lemma 4: For any vector /x = [hi,-- - , fi m ] with non-negative components, the optimal solution 
[S*, • ■ ■ , S£J for the following optimization problem: 

m 
i=l 



max 



subject to Ri = log 1 + 



V 1 ' :.,,;M,>W 
tr{Si)<Pi, SihO, i = l,---m, (31) 
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is also an optimal solution for the problem 

max h^Sihi 
subject to hjjSih 



13 ~ Hji 

tr(Si) < Pi, S t hO 
i,j = !,■■■, m, i^j, 



(32) 



*2 



with zfj = z tJ 



Following the same problem reformulation procedure used in Section II, to characterize the SUD rate 
region of an m-user MISO IC, the key appears to be the solution of (32) where zfj is a preselected 
constant denoting the interference power at the jth receiver caused by the ith transmitter, and Sj and Pi 
are the covariance matrix and power constraint for the ith transmitter. 

Unlike problem (4), the optimal covariance matrix for (32) with any given z^ 's need not necessarily 
be a beamforming matrix. Here is an example. 

Example 1: Consider the channels 

h xl = [1.9574,0.5045, 1.8645, -0.3398] T , 
h 12 = [-1.1398, -0.2111, 1.1902, -1.1162] T 
and h 13 = [0.6353, -0.6014, 0.5512, -1.0998] T . 

The optimal covariance matrix with the constraints P = 1, z\ 2 =0.3 and zj 3 = 0.6 is 

0.3784 0.1355 0.4103 -0.0461 
0.1355 0.0499 0.1587 -0.0220 
0.4103 0.1587 0.5443 -0.0965 
-0.0461-0.0220-0.0965 0.0273 
By restricting to beamforming, the optimal covariance matrix is 

0.2320 0.1198 0.3985 -0.0711 
0.1198 0.0619 0.2058 -0.0367 
0.3985 0.2058 0.6843 -0.1220 
l -0.0711-0.0367-0.1220 0.0218 

We have 



St = 



rank(S^) = 2. 



Si 



rank(Si 



1. 



h{S*hi = 7.1100 > h^li! = 7.0805 



t( 
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However, the above example does not mean that beamforming is not optimal for the SUD rate region of 
an m-user MISO IC. The reason is that the optimization problem (32) requires the interference powers 
to be exactly zfj. Without knowing the values of all z*- 2 , exhausting over all zfj will result in some 
rate pairs not on the boundaries of the SUD rate region. As we intend to establish the optimality of 
beamforming for achieving the boundary points of the rate region, we resort instead to the following 
more general formulation, i.e., the interference power is bounded by zfj, namely 

max hljSihu 
subject to hjjSihij < zfj, 

tr(Si) < Pi, Si y 

i,j = l,...,m, i^j, (33) 

This modified problem has its maximum that is no smaller than that of problem (32) for any values of zfj 
because of the relaxed constraint on the interference power. Therefore, if we can prove that beamforming 
is optimal for (33), then beamforming must be optimal for problem (31) even if it may not be optimal 
for (32). Such a strategy has been used in Theorem 1 where we let tpi vary in [0, | — 9j\ instead of the 
entire interval [0, |] because only in the specified interval does the useful signal power increase as the 
interference power increases. Based on the modified optimization problem (33), we obtain the following 
theorem. 

Theorem 3: For an m-user MISO IC, the boundary points of the SUD rate region can be achieved by 
restricting each transmitter to implement beamforming. 

The proof of Theorem 3 is achieved in the following steps. We first introduce Lemma 5 which allows 
us to solve problem (40). In the process, we establish that the rank of the entire covariance matrix can 
be set to the same as that of its submatrices. The final step is to show that the optimal submatrix need 
to have a rank no greater than 1 , established via an extension of Sylvester's Law of Inertia. 

We first introduce Lemma 5. 

Lemma 5: Let x and y be two complex vectors with dimensions t\ and ti respectively, and K be a 
(ti + £2) x {h + h) positive semi-definite Hermitian matrix with tr(K) < P. If 



K = 



Kn 14 

K21 K22 



(34) 
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and Kn is a preselected t\ x t\ positive semi-definite Hermitian matrix, then 

< (VxtKnx+HyHV^-MKn)) 2 , 

and the equality can be achieved by choosing K = K*, defined as follows. 
1) When x j Ki ,x / and / 0, we have 





t 




X 




X 




K 




y 




y 



K* = 



- — a/x T Kn 

. llyllVxtKnX 
2) When x j Ki ix = and \\y\\ / 0, we can set 



Kn 



K n xy 



/P-tr(Kn) 

Ij/HV^Khx 
P-tr(Kn) , 

— iu7rT2 — yy 
\\y\\ 



K* 



K 



ii 



yP-tr(K n , T 



where 



I/loKn 



In = 



\\v\\ 

P-tr(Kn) . 

— iu7h2 — 
l|y|| 



with 



and Kfj 





-l)xl 














Q, 



A 




Q 



being the eigenvalue decomposition of Kn, and A being a strictly positive diagonal matrix. 
3) When = 0, we can set 

K n 



K* 







Moreover, for all three cases, we have 

rank (K*) < max {rank (Kn) , 1} . 



(35) 



(36) 



(37) 



(38) 



(39) 



The proof is given in Appendix B. Here are some examples of Lemma 5: 

Example 2: Let x and y in Lemma 5 be two scalars, K^Obea2x2 matrix with tr(K) < 2. 
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If x = y = 1 and K\\ = 1, then from (36), we have a unique K* = , and rank(K*) = 1. 

• If x = 1, y = 1 and K\\ = 0, then from (37) we have a unique K* = , and rank(K*) = 1. 

2 

, and rank(K*) = 1. However, 



• If x = 0, y = 1 and K\\ = 1, then from (37) we have K* 

we can also choose K' = I2 X 2 and we have rank(K') = 2. Therefore, the optimal K in this case is 
not unique. 



1 1 
1 1 



1 




, and rank(K*) = 1. However, 



• If x = 1, y = 0, and K\\ = 1, then from (38), we have K* = 

we can also choose K' = 12x2 and rank(K') = 2. 
It is shown from the examples that only when x^~Knx / and / 0, or Kn = and \\y\\ 7^ 0, 
is the optimal K* unique. Otherwise we can choose other optimal K*s which may not satisfy (39). 
Lemma 5 is useful for the following optimization problem: 







t 






X 




x 


max 




K 






y 




y 



subject to /ij(Kn) = 0, i = 1, 
#(Kii)<0, j = l, 
tr(K) < P, KhO, 



,n, 
,m, 



(40) 



where /ij(-) and gj(-) are fixed functions. By Lemma 5, we can convert the above problem into 

2 





\v\Wp- 


tr(K n ) 


hi (K n ) = 0, 




• ,n, 


9j (Kn) < 0, 


i = i,-- 


■ ,m, 


tr(Kn) < P, 


k y 0. 





(41) 

Problems (40) and (41) have the same solution. Once the optimal Kn for problem (41) is obtained, one 
can construct the optimal K for problem (40) from (36), (37) and (38). As shown by Example 2, the 
choices of (37) and (38) may not be unique. One can choose K's that are different from (37) and (38) 
and still achieve the same maximum. 

With Lemma 5, we prove Theorem 3 as follows. 
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Proof: By symmetry, it suffices to show that for the mth user, the optimal covariance matrix S^ for 
the following optimization problem satisfies rank(S^) < 1: 

max hl nm S m h mm 
subject to h ] mj S m h mj <z 2 mj , j = !,-■■ ,m-l, 



<"mj _ ^mj > 

tr(S m ) < P m , S m y 0, 



(42) 



where all the /i m /s are t m x 1 vectors. 



We first show that problem (42) can be written as 

(y/tfSnh + V\\h mm \\ 2 - \\hf ■ y/ 'p - tr(Sn)) 



max 



subject to fJjSnhj < z^j, j = 1, ■ ■ ■ , m — 1 



tr(Sn) < 



Sn >- 0, 



(43) 



where h and all the /ij's, j = 1, • • • , m— 1, are m x 1 vectors, Sn is an m x m matrix, and m is defined 

as 



m = min {t m , m — 1} . 



(44) 



Obviously, when m = t rn < m — 1, problem (42) is exactly the same as problem (43) if we choose 
h = h mm , Sn = S m and hj = h m j. Thus, we need only show the equivalence of problems (42) and 
(43) when fh = m — 1 < t m . 



Let the SVD of h ml be 



h m l — Ui 



ll^mlll 
°(t m -l)xl 



and define 



s« = U t 1 S m U 1 
and h ™!' = U|h mi , j = !,-■■ ,m. 



(45) 
(46) 
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Substituting (45) and (46) into (42), we obtain 



max ^m^L^mm 



subject to h^sWh™ < z 2 mv j = 2, ■ ■ ■ , m - 1, 

t 

S (i) 



raj 
\\h m l\ 







(t m -i)xi 



ll^mlll 

0(t m -i) x i 



^ z mli 



tr(s«)<P m , S«^0. 



(47) 



Consider h^ 2 an d l et 



h (D _ 

n m2 — 



h 



(1) 
rn'2 



1 o 

u 2 



h 



(1) 

Hi2 



2,---,t„ 







(t m -2)xl 



(48) 



where ^h^ 2 ^ is a (t m — 1) x 1 vector consisting of the second to the last elements of h^ 2 - The 

SVDof(h2) ^ 



2,-,t„ 



(i) 

7/(2 



2,-" ,*r, 







(t m -2)xl 



where U2U2 = I(t m -i) x (t m -i)- Therefore 







t 






1 







1 








u 2 _ 







u 2 _ 



Define 



S (2) 



and h 



(2) 
mj 



1 

u 2 

1 

u 2 



S (l) 



4 ra xt„ 



1 

u 2 



*»£]■> J = ,m. 



(49) 



(50) 
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On substituting (48), (49) and (50) into (47), we have 



subject to hWs$h%] < z 2 mv j = 3, • ■ ■ , m - 1 



•mj 



\\h m i\\ 
°(t m -i)xi 

l m2 



1(2) 



H^ml || 
°(t m -l)xl 



< z, 



ml) 



m2 2- t 



S (2) 



°(t m -2)xl 

tr(s£>)<P m> S^O. 



t 



l m2 



l m2 j 



2,-,t„ 







(t m -2)xl 



(51) 



We note that the above transformation does not change the form of the previously modified constraints 
(see the third lines of problems (47) and (51)). Now we continue the above procedure up to h mjTn _i. In 
the jth transformation, we keep the first j — 1 elements of ■ and apply the SVD of the remaining 
(t m — j + 1) elements, and update the optimization problem. We formulate the jth, j = 2, • • • , m — 1, 
iteration as follows: 



.0-1) 



■JJ) 

mk 



h 



(j-i) 



I (i-i)x(j-i) 







0-i) 

"mj 



1,- J-l 



mJ / 7 ... t 

0(t m -j)xi 



L 0-l)x(j-l) 




u, 



and S# 



I (j-i)x(j-i) 







u, 



s 0-i) 



I(j-i)xO-i) o 

u, 



where (/^j denotes the jth to the i m th elements of 1 \ and its SVD is 



\ mj J 



= u,- 



v mj /,-,... 







where UjUj 



l (t m -j+l)x(t m -j+l)- 
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Finally, we convert problem (42) into the following form: 

h 
h 



max 



subject to 



h, 



- 1 


h 


s 






h 


1 1 





hi 



<z mj , j = l,---m-l, 



tr(s) 



< P m , sy o, 



(52) 



where h and hj are (m — 1) x 1 vectors and h is a (t m — m + 1) x 1 vector. Furthermore, \\h v 



\h\\ 2 + 



. Let 



Sn S 21 
S21 S22 

where Sn is an (m — 1) x (m — 1) matrix. The quadratic constraints in problem (52) are 



(53) 







= hJSii^ j <z^ i , j = l,---m-l. 



Therefore, the quadratic constraints in problem (52) are related only to Sn. By Lemma 5, problem (52) 
is equivalent to problem (43). 

Thus, in summary, problem (42) is equivalent to problem (43) with all the vectors in (43) being fhxl 
and Sn being fh x in. 

By Lemma 5, we can construct S in such a way that rank ^S^ < max {rank(Sn) ,l}. Let be 
optimal for problem (43). To prove Theorem 3, it is equivalent to prove 

rank(s^) < 1. (54) 

Furthermore, since S*i is optimal for problem (43), S* x also maximizes h^Snh under all the constraints 
in (43) with an extra constraint tr (Sn^ < tr (^11) • Therefore, to prove (54), it suffices to prove that 
the rank of the optimal covariance matrix for the following optimization problem is no greater than 1: 

max h^Snh 

subject to h^Snhj < z 2 m ^ j = 1, • ■ ■ , m — 1 

tr(Sn)<P, Sn^O, (55) 

where 



P = tr ( S* u ) < P. 



(56) 
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The equivalence is due to the fact that the optimal for problem (43) is also optimal for problem 
(55) and vice versa because of (56). Moreover, since is also optimal for problem (55), the inequality 
constraint tr (Sii^J < P is active. 
The Lagrangian of problem (55) is 



m— 1 



L = -rfSnh + X j (f^Suhj - z 2 m ^j + A m [tr (Sn) - P 



-tr (WSi; 



(57) 



j'=i 



On setting 



dL 



dS 



0, we have 



ii 



m— 1 



where 



W = -hh ] + Y \jhjh) + Ami 

3=1 

= C + X m I 



C = H*diag[-l,A 1 , - - ,A m _!]*Ht 



(58) 



(59) 



and H = [h, h±, ■ • ■ , h m -i] is an m x m matrix. 

We then introduce the following lemma which is an extension of Sylvester's Law of Inertia. 

Lemma 6: [29, Theorem 7] Let H be an s x t matrix and A be an t x t Hermitian matrix. Let tt(-) 
and v(-) denote, respectively, the numbers of positive and negative eigenvalues of a matrix argument. 
Then we have 

7r (HAH f ) < 7r (A) , and v (HAH*) < v (A) . 



By Lemma 6 and the Karush-Kuhn-Tucker (KKT) conditions that require Aj > 0, i = 1, 

v(C) < 1. 

Since C is an m x m Hermitian matrix, we can write its eigenvalue decomposition as 

C = Qdiag(r/i,-- - ,r/m)Q f 
where Q^Q = I, and r/j's are the eigenvalues in ascending order. From (60), we have 

m <o 

and rjj >0, j = 2, • • • , m. 



, m, we have 
(60) 

(61) 

(62) 
(63) 
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Since 



W = C + A m I = Qdiag [A m + r) U A m + 772, ■ ■ ■ ,X m + Vm] Q f , 

A m + rji is an eigenvalue of W, i = 1, ■ ■ ■ , to. 

Since the optimal S* x for problem (55) satisfies tr ( S* : ) = P, from the KKT conditions we have 



X m > 0. (64) 
Using (62)-(64) and noticing that the KKT conditions require W y 0, we have 

X m + m > (65) 
and A m + rjj > 0, j = 2, • • • , m. (66) 
For the optimal S^, from the KKT conditions, we have 
tr (wSJi) = tr ((C + A m I) S^) 

= tr (^Qdiag [A m + 771, A m + 7/2, ■ ■ ■ , A m + r/ m ] Q+Sn) 
= tr (^diag [A m + r?i, A m + rj 2 , ■ ■ ■ , X m + r/ m ] Q^^q) 

fa 

= ^(A m + r ?l ) (Q^JxQ).. 
= 0. 

Since (q^S^q) > 0, using (65) and (66) we have 

(Q t St 1 Q) ii >0, 

(Q^hQ) =0, i = 2,---,TO, 
v ' 33 

and ^Q^S^Q^ can be non-zero if X m + r\\ = 0. Since Q^S^Q >z 0, if one diagonal element is zero 
then all the elements on this row and this column must be zero. Thus, we have 

rank (S m ) < rank (S^) = rank (q^q) < 1. 



Theorem 3 proves the sufficiency of transmitter beamforming for achieving the SUD rate region. 
However, it does not mean that the SUD rate region can be achieved only by beamforming. As shown 



in the examples, the optimal S m is not unique when \\h\\ = or 
case in which h mm is orthogonal to h rn j, j = 1, ■ ■ ■ ,m — 1, and 
which h mm is linearly dependent of h m j, j = 1, ■ ■ ■ , m — 1. 



h 



h 



0. \\h\\ =0 corresponds to the 
= corresponds to the case in 
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B. SUD rate region of m-user MISO ICs 

Section III- A establishes the optimality of beamforming for the general m-user MISO IC with complex 
channels when each receiver is restricted to SUD. In this section, we first use a three-user real MISO IC 
as an example to show how we can obtain the SUD rate region by Theorem 3, and then generalize it to 
m-user complex MISO ICs. 

For a three-user MISO IC, the optimization problem (42) can be written as 

max hftSho 
subject to hJShi < z\ 

h 2 Sh 2 < Z 2 

tr(S) <P, S h 0. (67) 
We first reformulate this problem into (43). Let the SVD of h\ be 



fci = Ui 







where U^Ui = I. We then update ho, h\ h 2 and S as follows: 



h ( 2 L, =Vjh 2 

and = Uf SiUi, 

Ki) 



\\hi\ 


h 2 2 



(68) 
(69) 

(70) 
(71) 



where h 2 \ = (h^^j , and h 22 is the remaining part of h 2 ' . Let the SVD of h 22 be 

h 22 = \J 2 



h 22 
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and U^U2 = I. Then we update ho, h\ /i 2 and S as follows: 

T 

.(2) 



h\ 



(2) 



.(2) 



and S^ 
where h\\ is a 2 x 1 vector. Define 





1 














U ho 







u 2 _ 












T 






1 







\\hi 







u 2 _ 













T 






1 














u 2 




h>22 



1 

u 2 

: S^ = 



UfSiUi 



Sn s 21 
S 2 i s 22 



^01 

ho2 

\\hi\ 


fl22 



1 

u 2 



where Sn is a 2 x 2 matrix. Then we obtain the following optimization problem: 



max 



^•oiSii^oi + 11^02 1| J P ~ tt (^n) 



subject to 





T 






IIM 




IIM 




Sn 


< 














1 T 





^21 
^22 



'11 



< *2 



(72) 



(73) 



/l21 

^-22 

tr (Sn) < P, Sn b 0. 

The solution for problem (73) is complex. In the following we obtain the SUD rate region by using 
the fact that the rank of the optimal Sn is no greater than one without directly solving problem (73). 
Therefore, instead of exhausting over all feasible z\ and z\ and collect all corresponding Sn, we exhaust 
over all feasible Sn. We let 

sin 2 tpi cos tpi sin ip± sin tp2 

cos sin ipi sin tp2 cos 2 ip± sin 2 ip 2 



'li 



(74) 
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where ipi and Y> 2 can be any values in [0,7r]. Then we have that the optimal S for this Sn is 

S (h ,h 1 ,h2,P,i/)i,'ih) 

Sn 



Ui 



Ui 



1 

u 2 



1 

u 2 

1 

u 2 



: — - /, T . — =«02"oi & n 



\ho2\\ v ft 01 Snft i 



1^02 || V /loiSllhoi 

T^-h02H ()2 



P-trfSiOr fT 



1 

o 



when h 01 Snh 01 / 0, /i 2 7^ 



^02 



Oil f|T— fl S ll i ft 02 

II 02 1 1 

r 1T £5 p - tr <Sn)r £ 5 
ft02lo»ii I, e ii2 ft02ft c 



-T 
l 02 



1 




S n 




1 




uT 



when /» iSii/»oi = 0, 



*02 



7^0 



when 



02 



= 0. 



It can be shown that 



h-22 



h i = 



11*21 

ll*ol 



COS #12 

sin 6>i2 

cos #01 
sin #oi cos i 



where 



9ij = Z(hi,hj) , i,j = 0,1,2, 

x / COS #02 — COS #01 COS #12 



and # = Z (^02,^22) = < 



cos 



sin #01 sin #12 
tt/2, 



, #0l/0,7T, #1 2 /0,7T 

otherwise. 



# is the angle between the projections of ho and /i2 in foi's orthogonal subspace. The signal power and 
the associated interference powers are determined by Sn, and from (73) we have 

h^Sho = P \\ho\\ 2 (\ cos #01 sin^>i + sin #01 cos#cos ipi sin^l + sin #01 sin# |cos ipi cos ^2^ (75) 

hJSh! = P||/ii|| 2 sin 2 Vi (76) 

and /i^S^ = -P ||*2|| 2 (cos#i2 sin^i + sin#i2 cos^i sin^) 2 • (77) 

Therefore, we have the following theorem. 
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Theorem 4: The SUD rate region of a three-user MISO IC is 

T \ "1 ' 
/t n S(/tii, fei2, fel3, Pl,4>l,1p2)hll 

1 + h 21 S(h 22 , h 2 l, h 2 3, P 2 , 1p3, ^4)^21 + fyfi S (/I33 , /I31 ,h 32 , P3, 1p5, ^6)^31_ 

^22 S (^22, ^21, ^23, -P2, ^3, -04)^22 

1 + /lf 2 S(/ln, hi2, hi 3 ,Pi,1pi,l/)2)hi2 + ^32 S (^33, ^31, ^32, Pi, 1p5, ^6^32 

hl 3 S(h 3 3,h 3 i,h3 2 ,P3, 05, ^6)^33 

1 + hi 3 S(hu,hi 2 , hi3, Pi,ipi,ip 2 )hi3 + h 23 S(h 22 , h 21 , h 23 , P 2 , ip3, ip^h^ \ ^ 

(78) 

Furthermore, the boundary points of the rate region can be achieved by restricting each transmitter to 
implement beamforming. 

In the following, we discuss some special rate triples of the above region. For simplicity, we assume 
< Z (hij,hik) < §, where i, j, and k are integers ranging from 1 to 3. 

• ZF beamforming rate triple. 

When = 02 = 0, (75)-(77) become 

Hq Sh = P\\h \\ 2 sin 2 O i sin 2 0, 
/if Sfci = 
and h 2 Sh 2 = 0. 

Therefore, the rate triple in set (78) with ipi = 0, i = 1, ■ ■ • , 6 is the ZF beamforming rate triple. 
That is, the beamforming vector is chosen to be orthogonal to the channel vectors corresponding to 
the interference links associated with the same transmitter. 

• Single-user maximum rate surface. 

When = § - 0oi and Y> 2 = § - 9, (75)-(77) become 

/io S/i = P||^o|| 2 , 

/ifS/ii = P||^i|| 2 cos 2 O i 

and /ifS/^ = P||^2|| 2 cos 2 O 2- 

Therefore, h^Sho achieves the maximum with the constraint tr(S) < P. Then we obtain the 
following result. Let 

7T 7T _ : /cOsZ(hn,hi 3 ) - COsZ(/ln,^i 2 )cOsZ(/li2,/»13)\ 
V'l = 77 -Z(/tn,/li 2 ), V2 = --C0SM — — 7— , 

2 2 V sinZ(ftii,fti2)smZ(fti2,Ai 3 ) / 

4 ! 3 = tp5 = 0, ^4 e [0, 7r] and ^6 G [0, 7r] ; 



U 



^e[0,7r] 
i=l,- .6 



#1 < log 



#2 < 7J 



#3 < 2 l°g 
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then the corresponding rate triples form the surface of the three-dimensional (3-D) SUD rate region 
with R\ at its maximum. Similarly, the maximum rate surface for users 2 and 3 can be obtained. 

• Projection of the 3-D rate region to the two-dimensional (2-D) rate region. 

Suppose, say, user 1 is silent. We can recover the 2-D SUD rate region formed by users 2 and 3 
obtained from Theorem 1. That is, the surface of the 3-D SUD rate region corresponding to R\ = 
is precisely the 2-D SUD rate region of the two-user IC consisting of users 2 and 3. 

We can similarly obtain the SUD rate region of an m-user complex MISO IC with m > 3 by Theorem 
3. The only difference is how we generate matrix Sn. We generalize it as follows: 

S n = P 77 t, 

(7)1 = e3Ul1 sin^i 



and (7)^ = e 3Wz sin^i 



i-l 



k=l 

where tpi G [0, n], uii G [0, 2ir), and j denotes the imaginary unit. In the case of a real MISO IC, we can 
simply choose u>i = 0. The rate region and the corresponding beamforming vectors can be obtained in 
the same way as in Theorem 4, and so the determination of their quantities is omitted. 

IV. Numerical examples 

In this section, we present the SUD rate region of two-user or three-user real MISO ICs. For ease of 
visualization, we take the convex hull of all the regions. 

The achievable rate regions for the symmetric MISO ICs are shown in Figs. 3 and 4. Here, a symmetric 
MISO IC refers to one with \\hi\\ = \\h4W, \\h2W = \\h3W, 61 = 62, and Pi = P2. In Fig. 3, the rate region 
shrinks as 9 varies from | to 0, which corresponds to the MISO IC with the interference link and the 
direct link varying from being orthogonal to co-linear. The SUD rate region becomes smaller than that of 
frequency division multiplexing (FDM) when the direct link and interference link exhibit increasing linear 
dependence. When 6 = ^ , neither of the transmitters generates interference to the other user. Therefore, 
the IC reduces to two parallel single-user channels without interference, and the ZF beamforming rates 
achieve the maximum. When = 0, the two transmitters generate the worst interference and this MISO 
IC acts as a single antenna IC. Hence, the ZF beamforming rates are 0. 

In Fig. 4, the rate region also shrinks as the interference gain a increases. However, even if a — > 00 
the rate region will not be worse than the tetragon OMAN. Points M, N are the extreme points of the 
rate region on the axes. Point A denotes the ZF beamforming rates, which are determined only by Pi, P2 
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and 9. Therefore, as a increases, the rate region becomes close to the tetragon. Point A also shows the 
advantage of multiple antenna systems over single antenna systems, since the achievable rate region for 
the latter case reduces to the triangle defined by OMN when a — > oo. It can be shown that point A 
falls inside the FDM region only if 



sin^ < 



VI + 2P - 1 



(79) 



i.e., FDM outperforms SUD when the interference link subspace is close to the direct link subspace and 
interference gain is sufficiently large. Since 



U/T+2P- l U/T+2P- 1 
lim \ — = 0, and lim \ — = 1, 

P^oo V P P^O V P 

for all MISO ICs with 6>i ^ 0, § and 6 2 / 0, §, SUD achieves larger rate region than FDM for large 
power constraints, while the reverse is true for small power constraints and sufficiently large interference 
gains a 2 . 

Fig. 5 shows the SUD rate region of a two-user MISO IC under an interference power constraint. 
When such a constraint is small, this corresponding rate region is included in the FDM rate region. Since 
neither user generates interference in FDM, when the interference power is a concern in the system 
design, one can choose FDM instead of SUD when the interference is restricted. 

Fig. 6 shows the SUD rate region of a three-user MISO IC with the power constraints Pi = 1, P 2 = 1-5 
and P3 = 2. The channels are 

2.7 -0.5 1.2 1 

0.4 0.4 0.2 0.8 0.9 -1.7 

H, 1.5 0.9 ().;] ■ H 9 = -0.9 -1.3 -0.6 and H 3 = -2.6 0.8 -1 

0.8 0.4 0.3 1.3 0.7 

0.1 0.5 0.4 0.8 1.2 -1 

where Hi = [hu, h± 2 , hi 3 ], H 2 = [^21,^22,^23] and H 3 = [^31,^32,^33]- The solid curves are the rate 
regions for one user being inactive or at the maximum rate. That is, they are the projection of the 3-D 
rate region onto a 2-D plane with one rate fixed at a constant value. The ZF beamforming rate triple of 
this channel is shown in Fig. 7. The two cutting planes that pass through this rate point show that the 
ZF beamforming rate point is not on the boundary of the SUD rate region. 

V. Conclusion 

We have considered MISO ICs where each transmitter is limited to a Gaussian input and each receiver 
is limited to single-user detection. By exploiting the relation between the signal power at the intended 
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receiver and the interference power at the unintended receiver, we have derive a new method to obtain the 
SUD rate region for the MISO IC. We have shown that the original family of non-convex optimization 
problems is readily reduced to an equivalent family of convex optimization problems whose closed-form 
solutions are provided. As a consequence of restricting each receiver to implement single-user detection, 
transmitter beamforming is shown to be sufficient to achieve all boundary points of the SUD rate region. 




Qj i 1 i j 1 

0.2 0.4 0.6 0.8 

Ri (nat/Hz/s) 



Fig. 3. The achievable regions for the symmetric MISO Gaussian IC with \\hi\\ = \\h4W = 1, \\h2W = \\h3W = -j^,0i = 62 = 0, 
and Pi = P2 = 6. A1-A4 denote the respective ZF beamforming rates for 8 = ~ to 9 — 0. Also plotted for comparison is the 
FDM achievable rate region. 



Appendix 

A. Proof of Lemma 2 

Since it is straightforward to show (16)-(19) when h\ and I13 are linearly dependent, we only need to 
prove (12)-(15) under the condition that h\ and are linearly independent. 
Define 

r £n at 



a 



(80) 
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Fig. 4. The achievable regions for the symmetric MISO Gaussian IC with \\hi\\ — ||/i4 1 — 1, \\fl2 1 = \\h3W = u, 81 — 82 = f , 
and Pi = P2 = 6. A is the ZF beamforming rate point for all the choices of a. Also plotted for comparison is the FDM 
achievable rate region. 



where S\\ = [Sij * s a non-negative real scalar, a is a column vector, and A is a Hermitian matrix. 
Since Si y 0, we have 



S11A y aa 



(81) 



On Substituting (80) and (11) into (4), we have 

2 ^ 



h\S\h 



/in S 11 + h 11 fta + h[ 1 a J /3 + /3 J A(3 (82) 
IhfSn = zf, (83) 
where \\h3W 7^ since hs is linearly independent of h\. Therefore the optimization problem (4) is 



it. 



and hX^ihj, 
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equivalent to 



max 



.2 



\M 



in 



+ h u (3 j a + htaJp + 0* A/3 



subject to „/ 1 ,.- A >- aa^ , 

WM 2 ~ 

tr(A) < Pi 



Let the SVD of f3 be 



= 1^ 



11/311 




where UgUt = I. We further define 



A = UT AU^ 
a = uta, 



(84) 
(85) 

(86) 



(87) 
(88) 
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Fig. 6. The SUD rate regions of a three-user MISO IC. 



ill 



11 



and a.\ = • 



(89) 
(90) 



On substituting (83), (87) and (88) into (81), we have 



11*311- 

We first assume h\\ ^ 0. Then under the constraints in problem (85), we have 



(91) 



ll*3| 



(a) Zf 



111 



ll*3|| 2 



(b) z\ 



iim 2 



(c) 
< 



IIM 2 



'11 



'11 



'11 



o 



11/311 




+ 



11/311 




11/311 




+ h n \\p\\a 1 + h{ 1 \\(3\\a\ + \\f3\\ 2 A 



ii 



+ h n \\P\\a 1 + h\ 1 \\(3\\a\ + \\(3\\ 2 [P 1 



IIM 2 
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Fig. 7. The ZF beamforming rate triple Ri = 1.8118, R2 = 2.2998 and R3 = 2.3077, and the cutting planes that contain the 
ZF beamforming rate point and are parallel to R3-O-R2 and R3-O-R1, respectively. 



(<>) 

< 



< 



IIM 2 



IIM 2 

IIM 



/in 



'11 



'11 



+ 2||/3|| h n |ai| + ||/3|| 2 [P, 



IIM 2 



+ 2||/3|| 



'11 



IIM 2 



P< 



IIM 2 



+ II/3II 2 [Pi 



IIMI 2 



+ \\P\\\Pi 



11*3 



where (a) is from (86)-(88); (b) is from (89) and (90); and (c) is from the facts that A ^ and 



A n < tr (^AJ = tr (A) < Pi - 
and the equality of (c) holds if and only if 

z 2 

An = Pl ~W 



iiMr 



(92) 



(93) 



(94) 
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Since A y 0, (93) and (94) imply that 



7 2 



IIM 2 



(95) 



Inequality (d) is from the fact that 

hucii + ft-^aj = Re (jiuaij 
< 2 /in |2i| . 

The equality of (96) as well as the equality of (d) holds if and only if /inSi is real, i.e., 

t 



hucti = ( /inSi 
Inequality (e) holds because of (91) and (93): 



I Si I < 



An < 



II^H ||^3| 

The equality holds under the condition of (95) and 



'Pi 



(96) 



(97) 



(98) 



3i 



\h,\\ 2 



P 



z i 



IIM 2 







(99) 



where the choice of S ensures that (97) is satisfied. Thus, in summary, the equality of (92) holds if and 
only if (95) and (99) are satisfied. 

From (87), (88), (95) and (99), the optimal a and A are 

h\iZi 



a = XJpa 



and A 



-aa 



/in 
t _ 



'Pi 



hs\\ 
Pi 



JL 

11/311 



/3/3 f 
IIMV II/3II 2 ' 



(100) 
(101) 



Therefore, the optimal covariance matrix for problem (4) is 



Si = U 3 



huzi 



hi] 



h 



ii 



'Pi 



M\ 2 II/3II 



h 3 \ 
Pi 



'Pi 



\h 3 \\ 2 ||/3|| 



/3/3 j 



lfc.ll 2 ; ||/3|| 2 



U 



7i7i- 



(102) 
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The maximum of problem (4) is 

/ max 

Zl 




When /in = 0, problem (85) is solved by letting 



P 



/3/3t 



(103) 



In this case ct does not change the value of problem (85). It only needs to satisfy (81). In Lemma 2, we 
choose a to satisfy the equality of (81) to make the optimal covariance matrix be rank-1. Consequently, 
the optimal Si and the maximum of (4) are respectively (12) and (13) with k = 1. 



B. Proof of Lemma 5 

We first consider the special case of \\y\\ = 0. Obviously, (35) holds and we can choose K* as (38) 
so that (39) holds. One can also choose K22 / which still achieves the equality of (35) but violates 
(39). Therefore, in this case K* is not unique unless tr(Kn) = P. 

Another special case is Ku = and / 0. Since K ^ Owe have K21 = 0. To achieve the equality 
of (35), we choose K22 = PyyV\\y\\ 2 - Therefore, rank(K*) = 1 and (39) holds. 

Next, we need only to prove Lemma 5 when Ku / and ||y|| / 0. Let rank(Kn) = rank (A) = r. 
Then 



Ku 

K21 



x 



Q 
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K 2 i 
K22 
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K 2 i 

A2 1 





Q k 21 



K 



22 



I rxr 
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(a) 
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X~2 




y 




Xi 
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X 2 
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Xi 
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y 





K 21 Qt 






A 2 
I 



Q K 21 
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22 



Xi 
X-2 

y 



I rxr D 1 " 

D K 22 



Xi 

y 



Dt 



(c) 



Fi| 


\\y\\ 
o 

ll^ill 



llyll 
o 



L rxr 

D K 22 



Xi 



o u. 



Dt 



l rxr 

D K 22 



U X1 





o u„ 



Fll 



lll/ll 





_L rxr 



|Xi| 



lll/ll 





= ||xi|| + 2 ||xi|| 1 1 2/ 1 1 Re 

< ||2i|| 2 + 2 ||xi|| ||y| 

(d) 



(ujDU :Cl ) ii ] + ||y|| 2 (utK 22 U 2/ ) ii 
-||y|| 2 (utK2 2 U ?/ ) ii 



UtDU Xl 



< ||Xi|| 2 + 2 ||xi|| ||y||v / i 3 -tr(Kn) + Hi/11 2 (P - tr(Kn)) , 
( = } (VxtK n x + Hi/11 V^-t^Kn)) 2 
where in (a) we define 

Si A 2 n 

Qx, 



QK 21 . 



Xi 




Ai 





^2 







I 



and Xi is r x 1 and X2 is (ii — r) x 1. In (b) we define 



Dt 




A 



(104) 



(105) 



(106) 



(107) 





I 

The lower part of the matrix on the right-hand side of (107) must be the all-zero matrix, since the second 
matrix in the quadratic form (104) is positive semi-definite. In (c), we let the SVD of xi and y be 
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xi = U ; 



and y = XJ y 



\\ x i\ 


If I 




Since 



U^DtUy 



h 0, 



we have 



and 



uJdu Xi u] / k 22 u ?/ _ 

UtK 22 U y ) ^ < tr (uj,K 22 U,,) = tr (K 22 ) < P - tr (K n ) , 



UjDU Xl 



ii 



< J fu^K 22 U 2/ ) < v 7 ^ - tr (K 



ii 



ii • 



Therefore (d) holds with equality when 



UjDU Xl 



and UlK 22 XJ y 



y/P-tr(K u ) 

°(t 2 -l)x(r-l) 

P-tr(Kn) 

(«2-l)x(t 2 -l) 



From (106) and (108) we have 
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Therefore, (e) holds. Thus, in summary, the equality of (105) holds when (1 13) and (1 14) hold. Therefore, 
(34) is true and the optimal K* satisfies 



rank (K*) = rank 



rank 



I U^DtUj, 
U{DU X1 uJk 22 U,, 



y/P-tx(K U ) 

l rxr 



VP-tr(Kii) P-tr(Kn) 




= r 



= rank (Ku) . 

In the following we obtain K that achieves the equality of (35) 
1) When ajtK n x / 0, from (113) and (114) we have 

VP-tr(Kn) 



D = U„ 



and K22 = U 





P-tr(Kn) 




~ t _ VP-tr(Kn) 



ui 



Fill • \\y\ 



P-tr(Kn) f 
iu7h2 W ■ 

\\y\\ 



From (107) and (117), we have 

K 21 = [D 0] 
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y/P-tr(K n ) 



\ x i\\ \\y\ 



||y|| v/xTKnX 



(119) 



where the last equality is from (115). 

2) When x^KnX = 0, K22 is still (118). However, U Xl can be any unitary matrix since \\xi\\ = 0. 
Therefore, there are different choices of K21 that achieve the equality of (35). We choose U Xl = I for 
convenience. Then, from (107) and (117), we have 



_ VP-tr(Kn) T 
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